Chiral skyrmions in thin magnetic films: new objects for magnetic storage 

technologies? 
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Axisymmetric magnetic lines of nanometer sizes (chiral vortices or skyrmions) have been predicted 
to exist in a large group of noncentrosymmetric crystals more than two decades ago. Recently these 
magnetic textures have been directly observed in nanolayers of cubic helimagnets and monolayers of 
magnetic metals. We develop a micromagnetic theory of chiral skyrmions in thin magnetic layers for 
magnetic materials with intrinsic and induced chirality. Such particle-like and stable micromagnetic 
objects can exist in broad ranges of applied magnetic fields including zero field. Chiral skyrmions 
can be used as a new type of highly mobile nanoscale data carriers. 
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In magnetic materials with broken chiral symmetry 
the structural handedness induces chiral Dzyaloshinskii- 
Moriya (DM) couplings [1 which stabilize two- and three- 
dimensional localized structures with a fixed rotation 
sense and nanometer sizes [2j [3] . Originally, they have 
been described as chiral vortex-like configurations, but 
they are smooth and stable, topologically non-trivial 
magnetization configurations and, therefore, can be iden- 
tified as skyrmions in the micromagnetic limit with con- 
stant magnetization modulus, |M| = const . These 
skyrmions differ from other axisymmetric patterns in- 
duced by external dipole-dipole forces (bubble domains 
in nanolayers [4] and magnetic vortices in magnetic nan- 
odots [5 ) . Importantly, chiral skyrmions can also arise in 
nanolayers of magnetic metals where they are stabilized 
by surface/interface induced DM interactions [6 . In com- 
mon (centrosymmetric) magnetic crystals such solitonic 
states are radially unstable and collapse spontaneously 
under the influence of the applied field or anisotropy [2] . 
This singles out magnets with intrinsic and induced DM 
interactions into a particular class of magnetic materials 
where nanoscale magnetic solitons exist [3]. 

Recently, observations of such chiral skyrmions have 
been reported in nanolayers of noncentrosymmetric cu- 
bic ferromagnets (Fe,Co)Si and FeGe [7] and in mono- 
layers of Fe with a strong surface induced DM coupling 
[8]. This experimental break-through is not only an im- 
pressive demonstration of a unique phenomenon: static 
solitons and formation of solitonic mesophases in a chiral 
condensed matter system [3 j. These experiments also 
constitute new avenue for magnetic data storage and 
spintronics technologies. Chiral skyrmions, as magnetic 
inhomogeneities localized into spots of a few nanometer, 
can be freely created and manipulated, e.g., in extended 
layers of magnetically soft materials. These countable 
objects allow to fulfill a key task in creating versatile 
magnetic patterns at the nanoscale. 

This letter formulates the basic micromagnetic theory 
for thin ferromagnetic layers with chiral DM couplings 
and presents equilibrium solutions of isolated skyrmions 
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FIG. 1: (a) Axisymmetric isolated skyrmion with the core 
diameter 2R S in a thin magnetic layer of thickness L. (b) The 
magnetization profile along the diameter cross-section 0(r) 
and the perpendicular magnetization m z (r) indicate a strong 
localization of the skyrmion core; a (ro,#o) is the inflection 
point and the core radius R s is derived from Eq. (|6l). 



applicable to a broad range of material parameters. This 
is the first step towards a detailed calculation of the prop- 
erties and behavior of chiral isolated skyrmions in layer 
systems. 

As a model we consider a thin layer of a uniaxial ferro- 
magnet with intrinsic or induced Dzyaloshinskii-Moriya 
couplings. The micromagnetic energy density of the layer 
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A(gradM) 2 - M H - K(M • n) 2 + w d + w D (1) 



includes exchange energy with stiffness A, uniaxial 
anisotropy with constant K (n is a unity vector perpen- 
dicular to the layer surface), Zeeman, stray-field (w^) and 
DM (wd) energies [3j|9]. The chiral DM couplings are 
written as antisymmetric differential forms 
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that are known as Lifshitz invariants [lj. 

Introducing spherical coordinates for the magnetiza- 
tion M = M(sin 6 cos -0, sin 6 sin 0, cos 6 ) and cylindrical 
coordinates for the spatial variable r = (r cos (/?, r sin y?, z) 
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FIG. 2: Equilibrium magnetization profiles 6(r) for model 
([3| with if = 0, ft = 0.4, L/xq — 1 and different values of Q 
demonstrate a strong localization of the skyrmion core. Inset 
shows the reduced skyrmion core sizes R as a function of the 
reduced layer thickness L/xo derived by minimization of &(R) 
(Eq. @). 
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FIG. 3: The magnetic phase diagram in reduced variables k 
and applied magnetic field H/H a for fixed values of Q and the 
reduced layer thickness L/xo indicate the existence region of 
isolated skyrmions. In the double-hatched area spatially mod- 
ulated phases (helicoids and skyrmion lattices) correspond to 
the equilibrium state of the film. 



one can show that the equations minimizing functional 
(2) include axisymmetric localized solutions of the type 
ip(ip) (Fig. [TJ. The solutions ^(0) depend on the 
magnetic symmetry [2]. In this paper we consider stray- 
field free configurations ip = <p + 7r/2 arising in cubic 
helimagnets (crystal classes 23 (T) and 432 (O)) and in 
uniaxial ferromagnets with n22 (D n ) symmetry (n = 3, 
4, 6) [2J|3]. For these handed ferromagnets the part of 
the Dzyaloshinskii-Moriya energy responsible for the sta- 
bilization of skyrmions (i.e. that with in-plane gradients) 
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can be written as wjj = D(Ax. 
the Dzaloshinskii constant [2 j. The total energy of the 
axisymmetric string in the layer of thickness L can be 
written in the following reduced form W = 2ir ALw 
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sm 2 e + (H/H a )(l- cosO) +w d (p,L)/Q pdp. (3) [To] 



Here, we use a new spatial variable p = v/xq and char- 
acteristic parameters 
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where xq is the Bloch wall thickness, H a is the anisotropy 
field, Q is the quality factor [9], and the parameter n de- 
scribes the relative contribution of the DM energy. For 
k > 1 chiral modulations in form of helices or skyrmion 
lattices become equilibrium states in bulk magnets [2]). 
The stray-field energy of the axisymmetric string Wd is 
derived by solving the corresponding magnetostatic prob- 
lem [12 : w d (p,L) = (1 - 2sin 2 ((9/2)^(p,L)) 
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where S(p, x) = J °° Jo(&)Jo(pv) [1 - exp (-vx)} dv, 
Jq(x) are the zero-order Bessel functions. The equations 
minimizing the functional W ([3| with boundary condi- 
tions 0(0) = 7r, 0(oo) = yields the magnetization pro- 
files for skyrmions 0(r) (Fig. [I]) in the phase space of the 
four control parameters, k, H/H ai Q, L/xq. 

Typical solutions 9{r) derived by numerical minimiza- 
tion of energy functional ([3| are presented in Fig. [2] In 
a broad range of the control parameters the skyrmion 
profiles 6{r) consist of strongly localized arrow-like cores 
with linear variation of the angle ((tt — 6) ex r) and expo- 
nential "tails" with ex ^exp(— r/xo). A skyrmion core 
radius can be defined as 



R s =r -0 (dO/dr)v 1 . 



(6) 



where (r*o,0o) are the coordinates of the inflection point 
a and (d0/dr)o is the derivative in this point (Fig. [I]) 



The results in Fig. [2] demonstrate a variation of the 
skyrmion structure under the influence of magneto-dipole 
forces (parameters Q and L/xq). This allows to adjust 
the skyrmion structure and size by the variation of the 
layer thickness L or the value of the quality factor Q. 

In the phase diagram in reduced parameters H/H a 
and with fixed values of Q and L/xo (Fig. [3| we indi- 
cate the existence area of isolated skyrmions. Obviously 
skyrmions exist even in very high fields without collapse. 
At low fields the existence region of skyrmions is bound 
by several critical lines (Fig. [3|. They remain stable 
in zero and negative field for k < 1. At a certain criti- 
cal "bursting" field the skyrmion cores expand into 
a homogeneous state with magnetization parallel to the 
applied field. For n > 1 skyrmions either condense into 
lattices on the transition line H* or strip out into a helical 
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FIG. 4: Energy E of isolated bubbles and skyrmions as a 
function of their sizes (a) and the corresponding magnetiza- 
tion profiles 6(p)) (b) for a centrosymmetric magnet (D = 0) 
and that with finite Dzyaloshinskii-Moriya energy (D / 0). 



structure at a lower field H e \\. All these exceptional fea- 
tures of chiral skyrmions rely on the topological and en- 
ergetic stabilization of their core structure by chiral DM 
couplings. Therefore, chiral skyrmions are fundamen- 
tally different from cylindrical bubble domains [9 j, which 
are intrinsically unstable and only arise by the surface 
depolarization and the tension of ordinary domain walls 
as an effect of the shape of a magnetized body. 

In Fig. Hlwe demonstrate how axisymmetric solutions 
for model (pi) with D = (thin (blue) lines) transform 
into solutions with DM interactions (thick (red) lines). 
For D = only solutions for cylindrical domains (bub- 
bles) exist as metastable states in a certain range of pa- 
rameters H/H ai Q, L/xo [12 . Usually bubble profiles 
0(r) consist of an extended core with = tt separated by 
a thin domain wall from the surrounding homogeneously 
magnetized area with (9 = (Fig. [| b) [9j HI [12]. In 
Fig. [4] a energy (|3| plotted as a function of the bubble 
core size E(r) (blue line) indicates a metastable solution 
for r = R\. Under the influence of DM interactions the 
profile of the energy density is modified and includes so- 
lutions for skyrmions with fixed rotation sense and finite 
radius r = R S1 and solutions for isolated bubbles, which 
may have different rotation sense of the magnetization 
(r = R2 and r = R' 2 ) (red lines). The coexistence of 
skyrmion and bubble solutions occurs in a rather nar- 
row range of the material parameters. Outside this area 
bubbles are unstable and only skyrmions can arise in the 
film. 

In order to to elucidate the physical mechanisms for the 
stabilization of chiral skyrmions in thin magnetic layers, 
a simplified semi-quantitative discussion is useful. The 
profiles 6(r) in the skyrmion center are linear. The strong 
localization of its core allows to use a linear ansatz = 
7r(l - r/R), (0 < r < R), = 0, (r > R) as a suitable 
approximation for skyrmion solutions. With this trial 



function the equilibrium sizes of the skyrmion are derived 
by minimization of the function 

$ = A(H, Q)v 2 - 2k(x /L)v + v 3 G(v)/Q, v = R/L, (7) 

G{y) = 2 Jo 1 cos(7rr)rdr J* cos 2 «/2) S (r, f , v) £d£ with 
function 5(r, £,x) introduced in Eq. A(H,Q) = 

[1 + 8(1 - 4/7r 2 )(H/H a ) + (2/Q)(l - 2/tt 2 )}. The lin- 
ear term proportional to k is crucial for the existence 
of solutions with finite R. The "stiffness" coefficient A 
includes uniaxial anisotropy, stray-field and Zeeman en- 
ergy contributions. Finally, the self dipole-dipole energy 
of the skyrmion G(v) introduces a dependence of its size 
on the layer thickness L (Fig. [2] Inset). 

Recent experiments in nanolayers of magnetic mate- 
rials with intrinsic (cubic helimagnets Feo.5Coo.5Si and 
FeGe [7]) and with induced chirality (Fe/W bilayers) [8]) 
report observations of bound skyrmion states (lattices) 
(regions with ft > 1 in Fig. [3|. Isolated skyrmions 
with R s « 45 nm have been observed in a Feo.5Coo.5Si 
nanolayer (L = 20 nm) in the applied field larger than 
the critical field H* w 50 mT (Fig. [3| [7]. By using ex- 
perimental data from Refs. [3[T3] we found that for this 
sample n = 1.75 and xo = 16 nm. Our model results cor- 
roborate the identification of the observed magnetization 
patterns as chiral skyrmions. 

In conclusion, in magnetic layers with intrinsic or in- 
duced DM interactions isolated skyrmions with well- 
defined sized can exist as regular solutions of micro- 
magnetic equations in a broad range of the material pa- 
rameters (Figs. [2J [3]). Chiral skyrmions as particle-like 
spots of reverse magnetization can be considered as the 
smallest conceivable micromagnetic configuration. Im- 
portantly and in contrast to alternative and traditional 
storage technologies based on highly coercive or pat- 
terned media, chiral skyrmions can exist in magnetically 
soft nanolayers even at zero field, when sufficiently strong 
DM interactions can be induced them. In such materials, 
chiral skyrmions can be easily induced, transported, and 
controlled, e.g., by electric currents and applied magnetic 
fields. Stable chiral skyrmions in extended layer systems 
are promising objects for novel types of magnetic data 
storages, but also for logical bit- wise operations in ex- 
tended layer systems. 
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